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surface.  Next,  a  Gaussian  surface  neight  distribution  eer.eruUv  ■vsult  in 
less  diffuse  scattering  than  an  exponential  one.  The  trends  seer,  in  the  results 
for  Scattering  from  the  glistening  surface  do  not  differ  for  horizontal  or 
vertical  polarization,  for  Gaussian  or  exponential  surface  neight  distribution, 
or  for  different  signal  frequencies  in  the  S-band  to  I. -band  regions.  In  an 
absolute  sense,  the  results  for  vertical  polarization  indicate  that  there  would 
be  greater  angular  tracking  accuracy  tor  that  case  than  for  horizontal  polari¬ 
zation,  a  result  of  particular  interest  for  radar  applications. 
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The  Need  for  an  Expanded  Definition  of  Glistening  Surface 


.  IVIKOlUrilON 

Kleotrom.ignetie  signals  sc  <tt<  red  Irom  rough  terrain  include  contributions 
from  clutter  return  and  multipath  return.  ””  These  two  aspects  can  be  described 
thcoreticall’  it  properties  of  tne  terrain  such  as  the  probability  density  function 
li'lJ!  i  for  me  snrtaee  lieight  distribution,  tne  covariance  matrix  li,  the  variance  in 
-i: :  face  hctgnt,  "  and  the  complex  dielectric  constant  characterizing  the  surface  aie 
Known,  i’lf  uunit  runs  theoretical  models  of  KiM  wave  scattering  from  rough  sur¬ 
faces  * '  "*  all  relate  tile  norma li/.ed  cross  section  of  terrain  to  the  foregoing 

;  i :  ::,i  t.  r.->  >1  l : i •  ■  rough  surta.  e. 

II  UirkgM’lllnl 

Pile  problem  of  scattering  from  terrain  is  a  primary  interest  at  KAIX!  /  lsK( !. 
Taore  iiavi  be"n  number  of  models  developed  which  have  represented  the 
phenomena  involved  in  the  scattering  in  different  wavs.  Two  main  elements  common 
to  all  the  techniques  are  tile  description  of  terrain  features  by  the  use  of  statistical 
estimation  'iieorv  mi!  the  electromagnetic  scattering  formulation  that  incorporates 


i .  el .  e'l  a  . ,  ■  j  ubl  ii  a  lien  .1  anu a  r\  1  " : :  p 

!.),,■  a,  tile  iai'ae  number  ol  refe  renees  died  above,  tllev  .’.ill  not  he  listed  here. 


Extensive  descriptions  of  the  various  statistical  and  electromagnetic  models 
can  be  found  in  earlier  reporting.  ~  In  those  models  we  have  used  the  concept  of 
"glistening  surface"  (that  part  of  the  rough  surface  from  which  reflected  waves  can 
reach  the  receiver  for  a  given  position  of  the  transmitter  and  receiver'  as  des¬ 
cribed  by  Beckmann  and  Spizzichino.  ^  Their  expressions  for  the  length  and  width 
of  the  glistening  surface  were  derived  by  making  certain  assumptions  about  the 
physical  mechanisms  occurring  when  EM  waves  are  scattered  from  a  rough  surface. 
Some  general  discussion  of  the  models  will  be  outlined  here;  further  details  are 
presented  in  Appendix  A.  A  computer  program  has  been  developed  which  has  the 
capability  of  calculating  the  coherent  (specular  plus  direct)  and  incoherent  multi- 
path  power  reaching  a  monopulse  receiver  from  a  beacon  located  over  rough  terrain. 
Discussions  of  coherence  and  the  methods  for  calculating  the  coherent  and  incoherent 
power  are  given  in  Appendix  A.  The  results  that  were  obtained  are  discussed  in 
terms  of  comparisons  with  test  results  at  the  Discrete  Address  Beacon  System 

g 

(DABS'  site.  A  geometrical  representation  of  the  experimental  conditions  is 

shown  in  Figure  1. 
g 

Bahar  has  obtained  very  general  expressions  for  the  normalized  cross  sec¬ 
tions  (0°)  of  rough  surfaces  using  the  full  wave  approach.  He  has  shown  that  when 
the  major  contributions  to  the  scattered  fields  come  from  regions  of  the  rough 
surface  around  the  stationary  phase  points,  the  full  wave  solutions  reduce  to  the 
physical  optics  solutions.  In  our  models,  the  expressions  for  0°  are  based  upon 
the  physical  optics  approximations. 

When  the  surface  heights  are  normally  distributed,  the  expression  for  0°  used 
in  this  study  is  one  derived  by  Hagfors,  10  Barrick,  11  and  Semenov.  * ^  When  the 
surface  heights  are  described  more  accurately  by  a  bivariate  exponential,  then  o° 
is  given  by  an  expression  derived  in  Buck  et  al.  ^  Two  previous  studies1'^’  have 
indicated  that  the  suriace  heights  for  certain  terrain  are  better  described  by  multi¬ 
variate  probability  density  functions  having  bivariate  exponential  marginal  densities. 
This  result  comes  from  the  use  of  the  hypothesis  testing  procedure  (see  Appendix  A) 
developed  to  characterize  height  data. 

Another  important  aspect  is  the  inclusion  of  the  effect  of  local  shadowing. 

Sancer  s  '  shadowing  function  is  used  for  normally  distributed  surfaee  heights.  An 
extension  ot  Brown  s^  work  accounts  for  -shadowing  in  the  case  of  exponentially 
distributed  surface  Heights. 


t  Due  to  the  large  number  of  references  cited  above,  thev  will  not  hi-  listed  here. 
See  Ifete  rent  erf,  page  da.  t 
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'GLISTENING  SURFACE  IN  THE  PLANE  OF  THE  EARTH 

Figure  1.  Specular  and  Diffuse  Forward-Scatter  From  a  Rough 
Surface  and  a  Representation  of  the  Glistening  Surface 


1.2  llu-  Problem 

In  the  previous  studies,  the  use  of  the  Beckmann-Spizzichino  definition  in  the 
application  of  the  models  to  the  DABS  test  conditions  resulted  in  underprediction  of 
the  azimuthal  angle  uncertainty  u(  and  the  associated  diffuse  power,  particularly 
at  short  ranges.  As  a  result  ."e  were  led  to  an  examination  of  the  length  defini¬ 
tion  for  such  conditions.  As  an  alternative,  wo  chose  to  consider  the  entire  extent 
of  the  antenna  separation  as  the  surface  contributing  to  the  received  diffuse  power. 
The  defining  relations  for  the  two  surfaces  of  integration  are  the  following  (see 
Figure  D: 

First  there  is  the  length  as  given  by  Beckmann  and  Spi/.zichino.  '  Thev  deter¬ 
mine  the  length  bv  calculating  the  location  of  the  two  end  points  of  the  glistening 
surface.  1  or  a  homogeneous  surl  ice,  the  distant  -.s  from  the  end  points  to  the 

transmitter  and  receiver  1 1 .  ,  I  a  i  e  based  on  th<-  respective  heights  ( II  II.  I  and 

1  -  ■>  I  A 

a  roughness  criterion  uj/T1  wnore  a”  is  (he  surf  hc  height  variance  ind  T  is  the 

correlation  length.  1  hen. 


whc 


wn h  =  2  o i T  . 

H  o 

In  the  extended  length  definition,  the  entire  distance  along  the  surface  between  tne 
transmitter  and  receiver  is  included  and,  effectively,  =  1..,  =  U. 

For  both  cases,  Beckmann  and  Spizzichino' s  definition  of  width  (U  i  of  the 
glistening  surface  is  used: 


Here, 

D  =  total  ground  distance  from  transmitter  to  receiver, 

Xj  =  distance  from  transmitter  to  point  on  glistening  surface, 

Xj  =  distance  from  receiver  to  poii  t  on  glistening  surface. 

When  the  extended  length  was  included  in  the  model  it  was  found  that  better 
agreement  with  the  experimental  data  could  be  obtained  with  this  version.  The 
experimental  data  was  taken  by  Lincoln  Laboratory  personnel  at  the  Discrete 
Address  Beacon  System  (DABS'  test  site.  The  data  used  for  comparison,  boro- 
sight  azimuthal  pointing  direction  as  a  function  of  range,  lias  been  processed  to 
remove  the  contribution  associated  with  the  fact  that  the  mechanical  and  electrical 
boresight  directions  do  not  coincide.  The  bo  re  sight  pointing  direction  data  was 
divised  into  range  bins  and  a  local  mean-azimuthal-pointing-direction  was  estab¬ 
lished  for  each  bin.  Then,  the  standard  deviation  about  the  local  mean  was  calcu¬ 
lated.  This  standard  deviation  represents  the  azimuthal  angular  uncertainty,  u(J 
of  their  monopuise  receiving  antenna.  :  following  generally  accepted  usage,  this 
report  contains  three  similar  symbols  having  distinct  meanings,  t  are  should  be 
taken  to  distinguish  o,  a  ,  and  a,.  ' 

The  angular  uncertainty  increases  as  the  amount  of  diffuse  multipath  power 
entering  the  receiver  increases.  An  expression  relating  cj()  to  me  diffuse  power 
scattered  into  trie  receiver  is  given  in  Appendix  A.  It  includes  contributions  from 
noise  in  the  svateni. 

f  igure  2  shows  the  diffuse  power  for  the  two  theoretical  models  and  the  results 
based  on  the  experiment.  idle  terrain  data  base  allows  comparison  only  for  ranges 
less  than  60  km.  1’he  extended  definition  of  the  glistening  surface  provides  better 
overall  agreement  with  the  experimental  data,  particularly  for  ranges  K  less  than 
1 5  km.  For  those  range  s,  the  Beckmann  and  Spi/./ ieliiiu •  ,l< -fin i 1 1< >n  gives  'he  e’sai' 
•h-«t  the  glistering  surmee  do.  s  not  .-xist  ,r.  l  ‘hns  there  is  r.  r.  f  u 


significant  a  become  narrower  with  peaks  centered  near  the  specular  point  loca¬ 
tion.  When  the  region  over  which  a°  is  large  is  not  included  in  the  glistening  sur¬ 
face,  the  result  will  be  correspondingly  small  diffuse  power  contributions  and 
angular  errors.  A  further  point  is  that  even  for  sufficiently  large  ct“  values,  ther 
are  likely  to  be  o°  contributions  from  the  surface  areas  excluded  by  the  Beckmam 
Spizzichino  Lj  and  values. 

These  various  findings  provided  the  motivation  for  a  study  of  the  implications 
of  using  the  conventional  length  definition  and  consideration  of  when  it  would  not 
be  satisfactory. 


For  '.he  ;'iiT>o« 


this  studv  «e  selected  a  si  >!'  nominal  conditions  which 


wf"?  us'-.i  in  tile  a ;  •rstv  <>f  ases.  In  these  representative  eases  we  consider 
homogeneous  surfaces  where  me  heights  are  as.su, •  >  ed  t  >  have  a  bivariate  exponen¬ 
tial  dependence.  The  surfaces  are  considered  to  he  cultivated  terrain  with  complex 
dielectric  constant:  80  *  |  9.0.  For  consistency  with  the  nonuniform  data  base 
cases,  these  results  are  also  terminated  at  a  range  of  55  km. 

A  wide  range  of  results  are  presented.  The  first  aspect  is  the  consideration 
of  the  generality  of  the  effects.  Frequency,  polarization,  surface  height  distribu¬ 
tion,  ami  antenna  pattern  an  all  allowed  to  vary.  The  next  topic  is  a  parametric 
study  illustrating  'he  factors  controlling  the  nee  iyi  the  extended  length.  Then, 
some  specialized  conditions  are  considered.  Finally,  the  overall  question  is 
assessed  and  additional  topics  which  will  be  investigated  in  the  future  are  outlined. 

2.  GENERALITY  OF TIIK  RESULTS 

Most  of  the  cases  treated  m  >ur  main  study  involve  conditions  appropriate  to 
■he  DA  14S  experiment.  In  this  first  part  of  •  no  r<  port  we  wii.  sh<rt  that  although  the 
results  were  first  obtained  for  those  c  onditions,  thev  a><  not  limited  •<<  them  and 
-’filar  behavior  can  lie  found  tor  a  wide  nine'  >:  addi- i •  >t:al  conditions. 

In  ''-us  sec  tion,  the  rough  svrfa  ••  i -  .one  n  cl  to  i,  Imogen*-  ,.nd  'he  mean 
height  7  0.0.  The  height  if  tr  tsismitter  •  ;  !20h  on  :  die  height  . the  re  - 

■  eiver  is  100  rn.  unless  specific  1  otherwise. 

2.1  Slid, nr  lit" igli (  1  Hstrilxit ion  Hlnl‘ 

The  hypothesis  testing  pro.-i  lur<  s  describe  1  *n  \pncu  lix  A  wore  us  -d  to  deter- 
rrme  wl.,  (her  the  terrain  heigh.:-;  at  1  no  ItAHS  :<•  v  site  more  nearly  'if  either  a 
biv  iriat  Gaussian  or  exponential  iistribution  tun  dion.  1  he  result.-  of  mese  tests 
ini'  i  »  botti-r  fit  to  the  hiv-fiate  exponential  it  r  i  h'r  rn  fun .  :  i'  ••  .  ns  a  result, 
f.o-n  cf  't-...  invest. igat:  ins  on  r.«  ■•ffe.  t  of  exten  nng  'he  I'o.gth  of  ;h«  glistening 
-  n  r f  o  'c  :  n  this  repo  r't  ar*  -  'on  hi  '  i  ‘or  c  x  nonce* :  i !  9 s :  r: hut  ion s  in  -if  .  e  heights . 
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i  i  ■  >nci"  i  •!!'.  lirtf'nr  .fa  c  heights  m  I  tin-  He  .:.n  nr.  f .  1  Soi/zi- 
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Figure  7.  The  Effect  of  Assuming  Gaussian  Surface 
Heights  on  (Jp  Behavior  for  Both  Definitions,  Uniform 
Surface,  T  158.  114  m 


Another  set  of  results  for  a  is  presented  in  Figure  9  (horizontal  polarization) 
for  the  nonuniform  surface  (DABS  data  base)  and  various  antenna  height  conditions. 
The  equivalent  vertical  polarization  results  are  shown  in  Figure  10.  Although 
there  is  little  effect  when  both  antennas  are  close  to  the  surface,  there  is  some 
when  both  are  at  a  distance  above  the  surface.  The  most  significant  change  occurs 
when  the  two  antennas  are  at  different  heights.  The  extended  length  definition  for 
that  case  is  considerably  more  affected  than  is  the  Beckmann-Spizzichino  case. 

The  overall  horizontal  polarization  behavior  when  the  glistening  surface  is  extended 
is  analogous  to  the  case  for  vertical  polarization.  Note  that  for  all  the  antenna 
configurations  there  is  less  azimuthal  angle  error  for  vertical  polarization  than 
horizontal  polarization. 
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A  decrease  in  azimuthal  pointing  error  implies  better  tracking  accuracy,  and 
this  asocet  has  been  studied  further.  A  computer  program  that  simulates  ground 
based  radar  detection  and  tracking  of  multiple  targets  was  used  to  compare  both 
polarizations  in  terms  of  relative  tracking  performance.  In  the  simulation,  an  M  i  l 
radar  is  characterized  as  a  one-dimensional  rotating  phased  array  with  a  three - 
pulse  clutter  cancellation  system.  The  automatic  tracker  includes  ano-,1  filter  and 
has  six-target  maintenance  capability.  The  environmental  aspects  include  log- 
normally  distributed  ground  clutter.  Ravleigh  distributed  noise,  specular  multipath 
and  terrain  screening.  Several  track  initiation  and  maintenance  algorithms  an' 
employed.  Monte  Carlo  techniques  are  used  to  calculate  errors  in  range  ami  angle 
estimation  of  each  target's  position. 

1  he  results  of  these  simulations  show  that,  in  general,  the  MT1  1  adar  c.u,  more 
easily  initiate  and  maintain  track  on  multiple  targets  when  the  incident  held  i- 
vertically  polarized  than  when  it  is  horizontally  polarized.  This  is  found  to  li  d  1 
'or  various  surface  roughness  parameters  and  several  values  of  complex  dielec;  ric 
constant  characterizing  different  classes  of  terrain. 


Tlie  left  (land  pair  of  graphs  shows  both  glistening  surface  diffuse  power  results 

■>  2  > 

for  the  two  frequencies,  with  u"=  10  in'"  and  T-  500  in-.  The  right  hand  pair  of 

figures  shows  both  glistening  surface  diffuse  power  results  for  the  two  frequencies, 
■>  > 

w  ith  0~  -  10  m~  and  T  =  100  ni.  The  general  similarity  is  clear. 


2.4  Antenna  Pattern  Klfeetn 


In  Figure  12  we  show  the  diffuse  scattered  power  that  results  when  two  alter¬ 
native  elevation  plane  antenna  patterns  are  substituted  for  the  UAHS  configuration. 

The  terrain  is  assumed  to  be  a  uniform  surface  and  two  levels  of  roughness  are 

■>  2  '>  ’ 

considered.  T  500  m  with  u~=  10  m“  and  a“  100  ni“.  The  figure  shows  the  in¬ 
crease  in  diffuse  power  as  a  function  of  range  when  the  glistening  surface  length  is 
extended  for  both  an  isotropic  elevation  plane  antenna  pattern  and  a  (  x  j  pattern, 
it  is  clear  that  the  results  are  equivalent  for  either  pattern  The  corresponding  set 
of  results  for  the  DABS  antenna  pattern  can  be  seen  in  Figures  12  and  14.  In  those 
figures  the  right  hand  side  illustrations  are  the  diffuse  scattered  power  results  for 
the  conventional  and  the  extended  lengths  respectively.  Although  these  ranges  ter¬ 
minate  at  55  km,  it  is  evident  that  the  same  results  are  present  in  this  instance  as 
occurred  for  the  other  two  patterns. 
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Figure  12.  Antenna  Pattern  Fffects  for  Both  Definitions  With  T  500  m  and 
o^=  10  and  100  ni^ 
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!'<>  complete  the  discussion  for  T~  500  m,  we  want  to  consider  whether  there 

is  .ms'  significant  difference  in  the  effect  of  local  shadowing  when  the  alternative 

definitions  of  the  glistening  surface  length  are  used  in  the  model.  There  is  no 

2  2 

effect  on  diffuse  power  or  angular  error  when  o“  <  100  m“  for  either  definition. 

)  ’) 

for  e'  100  m",  there  was  a  slight  decrease  in  cr?  for  K  30  km.  This  can  be 

seen  bv  comparing  figure  15  with  the  corresponding  eases  in  f  igures  13  and  14. 

Tnis  is  consistent  with  the  expectation  that  the  effects  of  shadowing  would  become 
■> 

more  pronounced  as  u-  increases,  since  tins  corresponds  to  increasing  the  surface 
slop.  >.  \<>  significant  difference  in  shadowing  effects  can  be  attributed  to  the  use 

of  the  particular  glisti  ning  surface  definition. 

•'* Z i t.H . ccc.i . caacc 


l  igure  15.  ( 'ompa rison  of 

Shadowing  liffects  lor  the 
Two  Definitions,  T-  500  m, 
•>  > 
a~  100  m “ 


I  -  .i  ■'  .f  «;  ;-nd  diflusi  powei  with  range  for  vari  .ble  o"  values  has  been 
a  :  !>  ■  i  o.i  "ti  ul  a  •  or  v.  I.,  t  ion  length.  We  next  want  to  consider  tile 

.  ■  ■  t  .  -i'li  -  .t  s . nd  valm  ,  T  150.  Il  l  m.  This  value  was  selected 

!!■  ■  pr.  rts  ''  of  .  erf  dts  for  sotm  ratios  of  o  I  that  are  equivalent  to  the 

•>  '•  *  5 on  This  constant  r  ttio  comparison  will  be  discussed  in  more  de- 

i 1  ’  ■  ■  '  a  'M  4.  i'li'  present  discussion  is  lust  a  general  assessment  of  this 

f  ••  -s'. Its  to:  ly.th  length  definitions  mote  that  these  results  have 
o'!  l  .  o  •  bl  tefly  Iti  M-ctlon  2.  U. 

i  e.'  tne  ii'-okit.-inn  and  Spi/./.n  uino  definition,  f  igure  4  shows  the  typical  in- 

:  so  m  ■;  iti'i  jiff  use  power  with  increasing  The  shift  in  cutoff  from  58  km 

2  2  2  2  2  2  2 
1 5  kit.  .’.o',  o,  ,  or  .  -  '■  mceascs  from  J  1  m  to  o  10  m  .  for  a  -  100  ni  . 

i  ' :  ti  ill-,  •  ppe  •  r  >  . a !  ttige  the  :a 
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l  si*  "I  1  iu  t-xu-ndou  U*r4»t:i  <.it*i itiii u «n  ivsi.ns  m  i:, •.  1 •.  .  r  i  - •  n x ? : 

..rid  .i/imutiial  iinoTtiint.'  i-unuiart'd  n*  Uu-  1  ■;« *vk:i:.i an  .»no  unu  *.  :'t»r 
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u"  I  in”  .mu  o~  1  0  m~  tsi-e  1  igutv  51.  I  •>:•  those  \  alues,  the  .=  ;>:>  i n  n:  ..b.-vnee 

« •  t'  ingul  ir  uilccrtamtv  it  some  ranges  is  removed.  file  dnfe  fence  is  most  notice- 

able  tor  <J~  =  1  m".  There  was  no  effect  for  antenna  separations  less  than  55  km  m 

the  Heckmann-Spizzichino  results  while  full  integration  snows  that  actually  there  is 

>  > 

always  angular  error  present.  Of  particular  interest  is  the  case  for  u---  100  m-. 

1- or  that  case,  the  two  definitions  have  different  results  onlv  for  II  c  id  km.  Tor 
greater  distances,  the  lleckmann  and  Spizzichtno  definition  yields  glistening  sur¬ 
face  lengths  that  encompass  the  significant  o  regions  and  generate  results  equiva¬ 
lent  to  the  full  integration  eases. 

We  close  the  discussion  on  the  behavior  at  T-  158.  1  14  n.  bv  examining  the 

effects  ot  shadowing.  Comparison  of  figure  0  with  Figures  4  and  5  confirms  that 

again  the  effect  of  shadowing  is  the  same  for  either  definition  of  glistening  surface 

length.  For  both  definitions,  the  shadowing  results  in  a  decrease  in  un  onlv  for 
2  2 

a  =  100  ni  and  then  only  for  R  20  km. 

A. 2  Results  lor  Fixed  cr“  and  V  ariable  T 

In  this  section  the  behavior  of  multipath  effects  are  considered  for  fixed  values 

•> 

of  <J~  with  varving  T  witii  the  same  set  of  conditions  as  in  Section  .1.  1. 

Figure  16  shows  the  changes  in  cm  as  a  function  of  range  lor  the  lleckmann - 

Spizziehino  definition  when  cj"  -  10  m“  and  T  is  increased  from  1  m  to  500  m.  As  1 

increases,  bevond  I  -  100  m,  a  cutoff  in  <7j;  occurs  which  then  appears  at  successive 

Iv  larger  ranges.  As  T  increases  for  the  extended  length  definition  (see  Figure  17  i 

there  is  a  similar  tendency  towards  a  decrease  in  (J(  at  short  ranges  as  1  increases 

but  no  cutoff  occurs  at  any  T  value.  The  analysis  of  the  difference  in  results  is 

similar  to  that  introduced  in  Section  3.  1. 

The  preceding  discussion  has  been  for  the  case  where  there  is  no  shadowing 

2  2 

in  the  model.  Figure  18  introduces  the  effect  of  shadowing  for  cr“  =  10  m  and  a 
range  of  T  values  .  For  T  1  m  and  T  =  10  m  there  is  an  appreciate  shadowing 
effect.  As  T  increases,  the  shadowing  effect  lessen?  and  is  significant  mostly  at 
longer  ranges  until  for  T  -  500  m,  there  is  no  effect.  The  behavior  is  similar  for 
both  sets  of  glistening  surface  definition,  i  his  can  be  seen  bv  comparing  the  two 
T=  10  m  curves  of  Figure  18  with  their  respective  unshadowed  ones  in  Figures  16 
and  17. 
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Figure  16.  ati  Behavior  for  the  Case  of  BSL 
2  “  2 

Definition  a  =  10  m  ,  and  No  Shadowing 
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Computations  also  were  performed  for  the  case  of  Cr  -  1  m“,  and  T  varving 

from  10  m,  to  158.  114  m.  The  behavior  of  <J  and  diffuse  power  vs  range  is 

2  '  2 

analogous  to  the  case  where  o  was  fixed  at  10  m  and  T  is  increased  from  100  m 
to  500  m.  The  coherent  scattered  power,  though,  is  closer  to  flat  surface  conditions 
since  <7^  is  less.  Thus,  the  <J()  curves  have  the  corresponding  oscillatory  structure 
superimposed  on  the  general  variation.  The  general  comment  can  be  made  tiiat, 
for  fixed  ct“,  the  effect  of  shadowing  decreases  as  T  increases  since  this  corresponds 
to  a  decrease  in  surface  slope  (u/T).  Also,  the  shadowing  is  more  pronounced  at 
longer  ranges  which  represent  smaller  grazing  angles. 
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For  the  Beckmann -Spizzichino  length  ■i*-!  ii:m<>n.  the  giis-<  :d:ig  >  -  he  )- 

smaller  as  increases  and  the  (J(  curve  develons  a  i-utoti'.  Uiii-n  le'li  heights  -.re 
increased,  the  length  of  the  glistening  surface  approach-  zero  an  !  ih.er--  -  no 
diffuse  power  nor  o  ,  at  anv  range,  lor  ester,  to  i  lei.  -•!»  »  is--.  «1 10.  '  h''"  - 

power  and  d  decrease  as  either  II  ^  or  II  .  is  in-  t->-ase  i.  1  his  resiU"  ■■  a c 
attributable  la  a  decrease  in  the  width  ■*-  the  glistening  -uri'a.-c. 

Figure  20  shows  the  behavior  of  the  diffuse  power  for  the  two  sur  acc  tvp.  ; 
and  different  relative  antenna  heights.  1  he  extended  length  definition  i  -  used  an  • 
no  shadowing  is  present.  The  graphs  on  the  right  of  Figure  20  show  'he  unil'r.rr. 
and  nonuniform  surface  results  for  the  dual  antenna  height  condition.  The  results 
are  quite  similar.  Those  of  the  uniform  surface  are  slightlv  higher  for  all  ranges, 
with  the  greatest  difference  being  at  the  vicinity  of  the  dip  at  12  km.  I  his  dip  is  uue 
to  local  vanishing  of  the  width  of  the  glistening  surface.  I- or  the  low  height  case 
shown  on  the  left,  the  two  curves  are  again  similar  in  behavior,  will  'he  uniform 
surface  generating  slightlv  more  diffuse  power.  The  c  urves  no  longer-  show  anv  local 
minima.  Lowering  the  antenna  heights  increases  the  value  of  >j  p- ;  Jt  range-  less 
than  17  km.  with  a  more  rapid  fall  off  with  range  bevond  that  point.  This  is  tvu  nist 
due  to  changing  the  specular  point  since  the  uniform  surface  case  shows  similar 
behavior. 
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Figure  20,  The  Effect  of  Relative  Antenna  Heights 
on  Diffuse  Power  for  Both  Surface  Types,  El. 
Definition  and  No  Shadowing 
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Figure  10  is  .  1 1 r- : [ >  1  < * t . ■  I v  in.iiiii.Mius  m  Figure  1  *.  ex* -t ■  ot  that  the  terrain  charac¬ 
teristics  '■orrcspiin'  I  to  an  a.  tual  topographic  data  base  (see  Appendix  A)  instead  of 

a  um:orn.  region. 

An  examination  >f  Figure  10  shows  that  for  both  definitions  there  are  abrupt 
>seillations  in  plots  of  o  ,  vs  R  when  both  antennas  are  near  the  surface.  These 
variations  are  due  to  rav  blockade  (global  shadowing)  of  the  specularly  reflected 
rav  due  to  the  unevenness  of  the  terrain;  the  blocking  results  in  abrupt  variations 
in  the  coherent  sum  signal,  and  consequently  in  if  .  Consideration  of  the  case  where 
both  antennas  are  close  to  the  surface  in  Figure  19  might  suggest  that  there  should 
be  several  oscillations  in  a  as  R  varies  from  10  km  to  55  km.  However,  this  is 
not  the  ease,  because  here  the  mean  height  of  the  terrain,  Z  ~  60  m,  whereas  for 
the  uniform  region  (Figure  19)  7  0  m.  It  will  be  shown  insertion  4.3  that,  for 
Z  -  60  m,  there  are  ortlv  one  or  two  oscillations  in  P  qjj  between  10  km  and  55  km. 

Next  we  consider  the  dual  antenna  heights  case.  For  the  Beckmann-Spizzichino 
length  case,  the  cutoff  in  cj  t  as  H.j.  increases  (see  Figure  19)  occurs  for  the  non- 
uniform  region  results  also.  Similarly,  extending  the  length  of  the  glistening 
surface  also  results  in  a  noticeable  change  in  behavior  of  a  (  for  the  dual  heights 
ease,  binallv,  lor  both  antennas  at  a  large  distance  above  the  surface,  the  extended 
length  result  for  the  average  value  of  a  ,  vs  R  is  approximately  the  same  for  the 
nonuniform  surface  as  for  the  uniform  region.  This  does  not  hold  for  the  Beckmann  - 
Spizziehino  case  since  there  is  now  a  significant  extent  of  ranges  for  which  is 
detectable. 

4.  1.  3  TUB  BFFHCTS  OF  SHADOWING 

The  topic  here  is  the  effect  of  shadowing  for  the  various  heights  and  surface 
conditions  when  both  definitions  of  glistening  surface  are  used.  The  top  row  of 
graphs  in  Figure  2  1  show  the  results  for  with  shadowing  for  the  two  types  of 
surface  with  both  antennas  close  to  the  ground.  Since  both  definitions  lead  to 
similar  effects  for  that  case,  only  the  extended  length  results  are  shown.  Compari¬ 
son  of  these  two  graphs  with  their  corresponding  unshadowed  results  in  Figures  10 
and  19,  reveals  only  slight  decreases  due  to  the  shadowing  with  a  slightly  greater 
effect  for  the  nonunitorm  surface,  l  he  second  row  of  graphs  in  Figure  21  show 
the  result  for  o^  with  shadowing  for  the  nonuniform  surface  when  both  the  antennas 
are  at  a  large  distance  above  the  surface.  Both  definitions  have  been  used.  Com¬ 
parison  of  these  two  results  with  the  corresponding  unshadowed  cases  in  Figure  10, 
confirms  that  there  is  no  noticeable  contribution  from  shadowing  under  either 
definition  for  this  antenna  configuration. 
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relative  heights  of  the  antennas  are  changed  (DABS  data  basei.  In  both  equal  height 
cases,  the  specular  point  and  Fresnel  zones  tend  to  be  centered  between  the  two 
antennas  while  the  dual  height  <-ain»  results  if  spc  ular  regions  close  to  the  receiving 
intenna.  Thus,  the  substantially  different  behavior  hevona  .‘>0  Kit'  in  that  instance 
t  •'presents  contribution;;  to  the  coherent  power  from  region  "aving  a  iiflorent 
(smoother)  surface  •  liaracteri/atson  Tom  the  other  rases,  to  'enoral,  the  three 
curves  show  reiativelv  rough  surface  behavior.  1  he  sigtdn.  ant  difference  being 
that  for  the  low  height  case  there  are  some  i  av  blockage  •  cc’.s  present  and  some 
possible  multi:  nth  contributions  at  short  ranges.  This  is  -p. :  ui  h.v  Figure  2d. 

I  here,  the  7  ■  do  m  case  shows  tve  rt  behavior  for  a  •.>>;•  t-u-  tat  surface  with  the 
antennas  navin-r  iie  mean  antenna  nr.  l  I'd  of  the  nonam!-  .rrn  .  a  -<•.  I  tie  uniform  sur¬ 
face  removes  la  •  :av  i>l  mkage  ef  fee’s  f,r  l.-.w  i-.tc  nna  ne.g:  :  . 

The  rem.o  two  curves  o'  I-  ;gur”  id  .-..‘low  '  u  efre-  •  at  height  van  - 

•>  .  •/  •; 
anci>,  tot-  the  antennas  dost  to  ’.tic  urti.n:  tv.  sur'ai  e  w.*h  0  at.  For  o  “  10  m" 

the  coherent  power  shows  lar  'i  mr.piitu  i.  •  -••••  ..tations  for  H  .  »  am  tin  .  outcast. 


or  the  same  conditions,  the  dual  height, 
•a  illations  being  small  and  nri  m  ■  -• 


1  tide  of  '  in- 


can  be  explained  in  terms  of  the  fact  that  the  amplitude  of  the  coherent  power  is  a 

complicated  function  of  cj  sin  >'/A  where  A  is  signal  wavelength  and  sin  =  (II... -II  .  >''R 

2  1  ’ 
Then,  for  a  given  a  .  a  lower  transmitter  height  increases  the  oscillation  amplitude 

and  decreases  the  frequency  of  oscillation.  Thus,  onlv  at  large  K  values  does  the 

effective  surface  roughness  (a  sin  ill  A  )  become  small  enou  h  to  cause  oscillations 

in  the  dual  height  case,  while  the  roughness  is  always  sti  fieientlv  small  in  the  low' 

height  case.  Note  that,  bevond  the  R  75  km  point,  no  further  oscillations  are 

shown.  This  is  due  to  the  fact  that  the  large  amplitude  and  low  frequency  of  the 

oscillations  are  incompatible  with  anv  additional  plotting  within  the  scale  of  the 

■>  2 

figure.  In  contrast,  the  curve  for  a“  100  m  demonstrates  that  a  sufficiently  high 
variance  in  surface  height  preserves  rough  surface  behavior  even  for  the  low  height 
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To  .summarize  the  results  oi  these  studies,  we  start  with  discussions  of  the 
(a,  T)  aspects,  hirst,  there  is  the  behavior  of  o5  as  a  function  of  the  surface 
distance  from  the  receiver.  It  is  lound  that,  for  both  types  of  surface  height 
distribution,  graphs  of  u1  vs  range,  for  a  fixed  transmitter  location  exhibit  a  peak 
at  the  specular  point  for  small  <J/T  ratio  (see  Figure  3).  This  peak  becomes  less 
pronounced  as  the  ratio  (u/T)  becomes  larger. 

Next,  for  both  exponentially  and  Gaussian  distributed  surface  heights,  the 
behavior  of  the  diffuse  power  A-* ^ ) I i ' l-'  a  tuncf*orl  onlv  of  the  ratio  (o/T).  The 

coherent  power,  P  's  a  function  only  of  (J/A.  The  azimuthal  angle  pointing 
error,  (Jf)  is  a  function  of  both  and  Ppip'p”  ^ or  al*  altitudes,  the  effects  of 

extending  the  length  of  the  glistening  surface  on  the  behavior  of  Pjjjpp  and  vs  K 
becomes  more  pronounced  as  the  ratio  o/T  becomes  smaller.  For  small  o/T  and 
both  the  dual  height  and  high  height  antenna  conditions,  differences  in  Ppjpp>  of  as 
much  as  50  dB  have  been  observed.  For  these  cases,  the  Beckmann -Spizzichino 
definition  of  length  may  be  such  that  the  glistening  surface  does  not  include  those 
regions  which  contribute  appreciably  to  <J°  (see  Section  1  and  Figure  3).  In  contrast 
to  the  Pjjjpp  behavior,  as  the  ratio  o/T  becomes  smaller  the  effects  of  shadowing 
become  less  pronounced  for  the  two  surface  height  distributions. 

The  second  area  of  discussion  is  that  of  the  two  glistening  surface  length  defi¬ 
nitions.  Consider  the  effect  of  antenna  height.  For  the  Beckmann -Spizzichino 
length  case,  the  length  of  the  glistening  surface  becomes  smaller  as  either  the 
height  of  the  antenna  or  the  height  of  the  transmitter  is  increased.  Hence,  using 
the  extended  length  definitions  changes  the  behavior  of  Ppjpp  and  all  vs  R  for  either 
the  dual  or  high  height  conditions,  while  both  definitions  lead  to  equivalent  results 
for  the  hut  height  ease. 

Some  length  results  do  not  depend  on  relative  antenna  heights  or  surface  height 
distribution.  First,  the  effects  of  extending  the  length  of  the  glistening  surface  are 
more  pronounc  ed  at  short  range's  than  at  long  ranges.  Second,  extending  the  length 
) t  the  gliste  ning  surface  has  no  direct  influence  on  the  effects  of  shadowing.  Also, 
tor  both  length  definitions,  the  effects  of  shadowing  are  more  significant  at  longer 
ranges  (smaller  grazing  angles). 

Another  Tstinct  area  is  that  ot  coherent  power.  For  small  and  moderate  values 
if  a,  there  are  iscillations  m  vs  R  due  to  constructive  and  destructive  inter¬ 

ference  between  the  direct  ray  and  the  specularly  reflected  rav.  I  he  oscillations  in 
p  I  increase  in  amplitude  and  decrease  in  frequency  as  R  increases.  1  he  result 
•  I  this  is  that  corresponding  oscillations  appear  in  the  variation  of  d  vs  R  if  there 
is  sufficient  diffusely  scattered  power.  A  different,  coherent  power  phenomenon  is 
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Although  !ti<*  for  m  .  >'■  me  shadowing  function  for  tin-  '  .v  ■<  tv  i.un 

!i  ;;ril>u:  ion  is  liif.-rent,  thes"  studios  show  •,  hat  the  et'!.-  t.-.  »'  ,dowm.'  >r 
md  in  the  two  instances  arc  qualitati vclv  similar  i!'  a. I  other  1st  raim-te  rs,  arc 
u'|ual.  J  his  sin'ilariiv  is  also  lound  in  comparing  the  res  tills  of  th<-  two  length.'-"  )f 
glistening  surface. 

The  final  area  of  discussion  is  the  polarization  o'  the  signal.  1  he  effects  of 
extending  the  length  ot  the  glistening  surface  on  t'()[|.  ,  and  a  (  lot  horizontal 
nolarization  are  completely  analogous  to  those  tor  verti  al  polarization.  The  studii 
lead  to  tne  general  assertion  that  there  is  less  diffusely  s.  altered  rower  and  hence 
greap-r  angular  tracking  accuraev  when  the  incident  fi-l  I  is  vei  r.  eola>  d  zed. 

!  he  results  of  the  corroborative  analysis  of  this  aspect  allow  a  Tirtbci  coin  lu.;»on. 
Mm  the  meets  of  specular  multipath  for  vertical  polar! zation  result.  ;r.  better 
'rackiny  performance  of  an  M  !  I  radar  svsteni.  This  w  is  'our  :  to  he  true  at  h  ah 
i.-hand  and  S-lxuul  Iro'iuencies. 

We  have  discussed  the  ranges  of  terrain  statistics  an  I  systejr  parameters  lor 
which  the  usual  definition  of  glistening  surface  length  may  lead  to  -'ignifii  ant  under¬ 
estimation  of  the  diffuse  scattered  power  arid  azimuthal  pointing  error.  1  he  mag¬ 
nitude  of  the  effects  asso<  i  ale  cl  with  different  radar  or  t<  rram  -hara  ’(eristics  has 
also  been  assessed.  A  number  ot  additional  aspect:;  remain  to  be  resolved,  in  (he 
present  study  the  width  of  the  glistening  surfac*  was  dele ric-ined  bv  'he  Ueckmanri- 
Spizziehiito  relation.  Inasmuch  as  the  length  definition  has  proven  to  be  inade]ua 
m  mariv  cases  it.  may  also  be  true  that  removing  the  width  onstrain!  mav  alter-  the 
rt  suits  of  the  scattering.  This  can  be  explored  bv  in<  ludtng  an  azimuthal  variation 
in  in  the  calculations.  Another  aspect  that  will  be  investigated  is  tne  degree  of 
sensitivity  of  the  results  to  changes  in  ihc  dioleetrp  -  onstant  assigned  to  git -  u 
surface  features.  This  would  correspond,  for  example,  to  consideration  of  the 
o' feet  of  moisture  content  in  the  soil. 
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Appendix  A 

Scattering  Model 

The  results  of  this  study  were  obtained  from  a  sophisticated  model  describing 
the  multipath  effects  associated  with  scattering  of  a  signal  by  rough  terrain.  This 
model  includes  such  factors  as:  la)  azimuthal  and  elevation  power  patterns  of  a 
monopulse  receiver  (see  Figure  Al);  (b)  the  spatial  nonuniformity  of  the  rough 
earth;  (c)  nonuniformities  in  the  glistening  surface;  (d)  possible  multiple  specular 
reflection  points  due  to  unevenness  in  the  surface;  (e)  the  surface  height  distribu¬ 
tion  characterization  and  (f)  both  global  and  local  shadowing  effects. 

The  program  normally  uses  system  parameters  associated  with  the  I. -band 
Discrete  Address  Beacon  System  to  allow  comparisons  with  experiment.  These 
are  defined  in  Table  Al  and  Figure  Al.  External  inputs  include  the  complex  dielec¬ 
tric  constant  for  each  surface  area,  the  coordinates  of  the  monopulse  receiver, 
the  velocity  and  initial  and  final  position  of  the  aircraft  containing  the  transmitter, 
and  a  parameter  to  control  the  effects  of  shadowing. 

The  description  of  the  model  may  be  divided  into  two  major  topics.  These  are 
the  techniques  required  to  assign  appropriate  statistical  properties  to  the  terrain 
and  the  specifics  of  the  electromagnetic  formulation.  Each  of  these  aspects  re¬ 
quires  some  discussion. 

There  are  several  surface  feature  contributions  in  F  model.  Analyses  of  the 
scattering  from  rough  surfaces  consider  the  surface  heights  in  the  region  as  pairs 
of  scattering  elements  and  in  most  cases,  including  this  study,  assume  that  the 
height  distribution  can  be  described  bv  either  a  bivariate  Gaussian  or  exponential 
probability  density.  These  two  bivariate  densities  have  the  forms: 
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simple  alternative  hvpothesis  test  is  equivalent  to  a  nitnimum  error  probability 
criterion  where  it  is  equally  likely  that  either  density  is  the  appropriate  on*'. 

I  or  the  comparison  of  tile  scattering  theory  with  the  experiment,  these  tech¬ 
niques  had  to  be  applied  to  the  terrain  at  the  eastern  Massachusetts  site.  A  data 
base  of  topographic  derations  for  this  area  is  available  at  the  Flectromagnetic 
Compatibility  Analysis  Center  « Ft  AO  in  Annapolis,  Maryland.  This  was  prepared 
from  digitized  terrain  maps  supplied  bv  the  Defense  Mapping  Agency  (DMA'.  The 
area  of  interest  is  divided  into  rectangular  cells,  each  with  sides  of  about  -  km. 
Much  cell  is  further  subdivided  into  a  10  bv  10  grid  of  points.  The  statistical  analy¬ 
sis  is  then  applied  to  the  individual  cells. 

The  statistical  data  for  each  cell  has  been  recorded  on  a  computer  tape  for  use 
with  the  program  for  the  electromagnetic  analysis.  Mach  cel!  is  represented  by 
seven  descriptors.  1'he  first  two  entries  are  the  (x,y>  coordinates  for  the  center 
of  the  cell.  Next  is  the  geological  code  t dielectric  constant)  for  the  cell.  (The 
predominant  feature  is  woods;  there  are  a  number  of  cells  containing  clusters  of 
lakes  and  ponds  and  a  few  town  sites  with  associated  cleared  areas.  )  This  is 
followed  by  the  mean  and  variance  of  the  heights  in  the  cell  and  the  estimated  corre¬ 
lation  length,  T  (the  units  of  length  are  in  meters).  The  final  quantity  is  the  result 
of  the  hypothesis  test. 

The  trajectory  of  the  beacon  aircraft  is  incorporated  into  the  computer  program 
and  at  each  range  point  for  which  a  calculation  is  to  be  made  the  required  cells  and 
their  descriptors  are  then  identified.  These  results,  or  in  the  general  case,  the 
equivalent  set  of  input  parameters  are  then  used  in  the  electromagnetic  analysis. 

The  calculation  of  the  electromagnetic  fields  has  two  distinct  elements.  First, 
the  total  coherent  electric  field  F  .  at  the  receiver  is  calculated  using  the  sum 
pattern  of  the  receiver  antenna  in  the  following  expression: 


y'fa.ifl  >K.  I"  elkAH] 

A-<  i  m  l  I-  ' 

.1  J  .1 


T 

G  AO  ) 
I  ni 


-  direct  path  electric  field  at  the  receiver, 

-  gain  of  receiver  in  direction  of  multipath  ray,  where  tt  ^  is 
the  angle  between  the  direct  ray  and  multipath  ray, 

-  attenuation  factor  affecting  coherent  reflected  wave  due  to 
surface  roughness, 

-  complex  Fresnel  plane  wave  reflection  coefficient, 

=  2  t;  /  A  ,  and 

-  difference  in  path  length  between  direct  ray  and  reflected  ray 


The  summation  over  j  represents  all  possible  specular  reflections  (there  may  be 
more  than  one,  due  to  unevenness  in  terrain).  Here,  coherence  means  there  is  a 
known  phase  relationship  between  the  Hl\. I  field  leaving  the  transmitter  and  that 
reaching  the  receiver. 

The  next  aspect  is  that  of  the  diffuse  power.  The  diffuse  power  in  the  mono¬ 
pulse  difference  channel  is  calculated  from  the  expression 
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system  processing  losses, 
transmitted  power, 
wavelength, 

gain  (power)  of  transmitter  in  azimuth  (isotropic  pattern), 

gain  of  receiver  in  azimuth  (difference  pattern.  Figure  Al), 

gain  of  transmitter  in  elevation  (isotropic  pattern), 

gain  of  receiver  in  elevation  (Figure  Al), 

elevation  angle  between  bo  resight  and  point  on  glistening 
surface  for  transmitter, 

elevation  angle  between  boresight  and  point  on  glistening 

surface  for  receiver, 

range  between  transmitter  and  point;  on  glistening  surface, 

range  between  receiver  and  point  on  glistening  surface, 

element  of  area  of  glistening  surface  which  is 
illuminated  bv  beacon, 

azimuthal  angle  between  boresight  and  point  on  glistening 
surface  for  receiver,  and 


azimuthal  angle  between  boresight  and  point  on  glistening 
surface  for  transmitter. 
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with  a  similar  expression  for  the  horizontal -horizontal  scattering  matrix  elemi 
fi  I  .  In  the  above  expressions,  the  subscript  .  ,  refers  to  the  I, -field  in  the  ;>!• 
of  incidence: 

(  cos  a  -  V'  -  sin”  « 

K  *nr  1  =.  _i - l—L - 

t.  (  cos  Of  -  y  i  -  sin^  ot 

<.  -  the  relative  complex  dielectric  constant  of  the  surface, 

t>  -  angle  of  incidence  (with  respect  to  surface  normal), 
d  .  ~  angle  of  scattering  (with  respect  to  surface  normal). 


tor  an  exponential  surface  height  probability  density  function.  The 
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